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In the context of an effective model for doped antiferromagnets, whereby the charge carriers are 
treated as hard-core bosons, we demonstrate that the ground state energy close to half-filling is an 
even periodic function of the external magnetic flux threading the square lattice in an Aharonov- 
Bohm geometry. The period is equal to the flux quantum $o = 2nhc/q entering the Peierls phase 
factor of the hopping matrix elements. Thus flux quantization and a concomitant finite value of 
superfluid weight Da occur along with metallic antiferromagnetism. We argue that the charge q in 
the associated flux quantum might be set equal to 2e. The superconducting transition temperature 
Tc is related to Ds linearly, in accordance to the generic Kosterlitz-Thouless type of transition in a 
two-dimensional system, signalling the coherence of the phase fluctuations of the condensate. The 
calculated dependence of Tc on hole concentration is qualitatively similar to that observed in the 
high-temperature superconducting cuprates. 
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I. INTRODUCTION 

The continuous evolution of the charge and spin dy- 
namics upon doping with mobile holes the insulating an- 
tiferromagnetic copper oxide layers, supports the view 
that the high-temperature superconductivity in these 
materials |l| is a fundamental property of the two- 
dimensional doped antiferromagnets 0. Of particular 
importance in this respect is the observed linear increase 
of the superfluid weight Ds with small hole concentra- 
tion (1 — Ue), away from the Mott metal-insulator tran- 
sition point at half- filling |g[. On the theoretical side of 
this problem, there has been some evidence that mobile 
holes in doped antiferromagnets behave much like hard- 
core bosons. This transmutation of statistics, from bare 
fermionic holes to bosonic vacancy quasiparticles, should 
be understood as an "emergent phenomenon" due to the 
reduced dimensionality and the presence of a strongly 
correlated spin background. In the context of the simple 
fermionic t-J model, proposed by Anderson [Q to de- 
scribe such systems, the aforementioned evidence comes 
from exact-diagonalization studies of the ground-state 
energy and the static hole-hole correlation function on 
small clusters (|-01- Indeed, the possibility of a hard- 
core boson behavior of the charged vacancies in doped 
antiferromagnets, opening the way to Bose-Einstein con- 
densation and the appearance of superconductivity, was 
suggested by many authors ||,y in the early days of high- 
Tc superconductivity research. Thouless ^, in particu- 
lar, argued that due to topological constraints, a vacancy 
in a two-dimensional torus lattice threaded by an exter- 
nal magnetic flux must be transported twice around the 
ring in order to recover its original configuration. Hence 
flux quantization with an effective charge q = 2e may 



result from this period-doubling of the charge e bosons. 

In all the aforementioned works, the lack of an effec- 
tive model for doped antiferromagnets expressed in terms 
of hard-core bosons has prevented the systematic study 
of their flux quantization properties in conjunction with 
the optical and magnetic ones. Such a model, however, 
has been postulated from the outset by Psaltakis and Pa- 
panicolaou iQ and consists of a t-t'-J Hamiltonian and a 
suitable 1/A^ expansion that provide a reasonably simple 
many-body calculational framework for the study of the 
relevant issues. When leading quantum-fluctuation ef- 
fects are taken into account in the context of this model, 
the generic experimental features of the optical conduc- 
tivity, the Drude weight and the total optical weight in 
the cuprates are qualitatively reproduced. In particular, 
our theory |l^,0] accounts aptly for the experimentally 
observed 0.5 eV peak of the midinfrared band |13|-p5[ and 
the mass enhancement factor approximately equal to 2 
p4| . Furthermore, it predicts a finite limiting value for 
the optical conductivity a{u! — > 0), at finite hole dop- 
ing, consistent with the residual far-infrared conductiv- 
ity observed in the YBa2Cu306+a; family of cuprates p6[ . 
Our results are also found to be consistent with relevant 
exact-diagonalization data Jl^l ■ 

In view of the quoted evidence from optical experi- 
ments in favor of our effective model, we have recently 



undertaken a systematic study |18| of its flux quanti- 
zation properties in order to provide a more complete 
assessment of the main electromagnetic responses. Our 
study includes results for the superfluid weight Ds and 
the associated superconducting transition temperature 
Tc- In particular, our explicit numerical estimates for 
the doping dependence of Tc, including leading quantum- 
fluctuation effects, are found to reproduce qualitatively 
the observed trends in the cuprates ||,|l^ . In the follow- 



ing we review the main points of this approach. 

II. EFFECTIVE MODEL 

Our effective model is described by a t-t'-J Hamil- 
tonian expressed in terms of Hubbard operators x"'' — 
\a){h\ as 












xt'xT)^ (1) 



where the index corresponds to a hole, the Greek in- 
dices ^,v, . . . assume two distinct values, for a spin-up 
and a spin-down electron, and the summation conven- 
tion is invoked. Here J is the antiferromagnetic spin- 
exchange interaction between nearest-neighbor sites {i,j) 
on a square lattice endorsed with periodic boundary con- 
ditions and a total number of sites A = A^; x Aj^, where 
A J, = Aj^. For the hopping matrix elements i^ we as- 
sume non-zero values, t and — t', only between nearest- 
and next-nearest-neighbor sites, respectively, as dictated 
by quantum-chemistry calculations i2Q;El| for Cu-0 clus- 
ters and fits of the shape of the Fermi surface observed 
by angle-resolved photoemission spectroscopy |22|. We 
also generalize the local constraint associated with P) to 
X™ + xt^ — ^j where N is an arbitrary integer, and 
consider the commutation properties of the x'''"s to be 
those of the generators of the U(3) algebra. A represen- 
tation of the latter algebra in terms of Bose operators 
can then be employed, leading to a generalized Holstein- 
Primakoff realization that resolves explicitly the local 
constraint which gives rise to the hard-core character of 
these bosons. One can then develop a perturbation the- 
ory based on the 1/N expansion, restoring the relevant 
physical value A^ = 1 at the end of the calculation. 

In the presence of an external magnetic flux <&, thread- 
ing the two-dimensional lattice in an Aharonov-Bohm 
torus geometry, the hopping matrix elements i^ are mod- 
ified by the well-known Peierls phase factor and should 
be substituted in M according to 



^ij ^^ f^ij 6 



with Aij ^ 



(R, - R,) • e, . (2) 



Here R.^ is the position vector for site i, e^ is the unit 
vector along the x-axis encircling the flux lines and 
$0 = '2,nhc/q is the so-called flux quantum pJ. Con- 
ventionally, the charge q of the carriers entering $0 is, of 
course, equal to the electronic charge e. However, the ar- 
guments of Thouless B quoted in the Introduction imply 
that a vacancy actually "feels" twice as much external 
flux. In the context of the present effective model this 
may be accounted for by an extra factor of two in the 
expression (0) for the Aij which can be readily absorbed 
in a redefinition oi q as q = 2e. Evidently, this reason- 
ing does not constitute a rigorous justification for the 



assignment q — 2e in the flux quantum $o- The latter 
justification can be provided only by an ab initio deriva- 
tion of an effective Hamiltonian for the hard-core boson 
vacancies, starting from a realistic electronic model for 
the cuprates. At present such a program is out of reach. 
Hence this work will be content with the study of the flux 
quantization properties of the effective model described 
by (^-(||), given the flux quantum constant <I>o. 

In the large- A'^ limit the Bose operators become classi- 
cal commuting fields which can be parametrized by the 
local electronic density ni, the angles 9i and 4>i deter- 
mining the local spin direction, and the local phase ipi 
of the condensate. The Hamiltonian (|^)-(0) takes the 
form H{^) = N'^AEoi^), where Eo{^) is the classical 
energy per lattice site for the value of physical inter- 
est A^ = 1. More explicitly, for uniform density states, 
Ui = He, where rie is the average electronic density, we 
have that 



AEoi'^) = £1 + £2 



(3) 



where 



£1^- 


-ne{l-ne)y^t,j 
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(4) 


£2^- 


-j-jy^ [cos 6i cos 6j + sin 9i sin 0j cos(0i — cfij ) 


-!]• 



As shown in Ref . [|T^ , close to half-filling (rig < 1) and for 
a sufficiently large t' , the ground state of (||)-(|j), in the 
absence of an external magnetic fiux (<I> = 0), is described 
by a planar spin configuration (6i — tt/2) in which the 
local twist angles and phases are modulated according to 



= Q R, , ^, = Q' R, 



(5) 



where Q = (tt, tt) is the usual spin-modulating anti- 
ferromagnetic wavevector and Q' = (tt, — tt) is an un- 
usual phase-modulating wavevector. The question that 
is now posed is how this phase-modulated antiferromag- 
netic (AF) ground state will respond to the presence of 
an external magnetic flux $. 



III. FLUX QUANTIZATION AND SUPERFLUID 
WEIGHT 



Following an argument by Yang [E4J we note that, in 
the presence of $, the reciprocal lattice is displaced from 
the origin by 2tt^/{Ax^o) along the x-axis. The quanti- 
zation of flux therefore depends on whether the ground- 
state energy of the system changes under this momentum 



boost. Given that the spin-exchange part of the Hamil- 
tonian (0)-(||) does not couple directly to the magnetic 
flux it is plausible that, at least in the large- TV limit, 
the condensate will respond in such a way as to leave its 
spin-modulating wavevector Q intact and simply adjust 
its phase-modulating wavevector Q' to a new value. In 
other words, we anticipate that in this classical (large- 
iV) limit, the rigidity of the ground state against the 
intrusion of the external magnetic flux comes solely from 
the phase fluctuations of the condensate. These heuris- 
tic arguments lead us to consider the ansatz (H) with 
the following modulating wavevectors: Q = (tTjTt) and 
Q' = (tt, — tt) — (47r?7i/Aa;, 0), where m is an arbitrary 
integer. Inserting these wavevestors into (|5|) and taking 
carefully the infinite lattice limit (A -^ cxd) of the classical 
Hamiltonian (§)-(Q) we have that 



KEq{<^) - AEq{<^ = 0) = 8i'7r^ne(l - Ue) 



$0 



(6) 



Thus for each integer m we get an individual many- 
body energy level that depends quadratically on <I>. The 
ground-state energy is given by the lower envelope of 
these crossing energy-level parabolas and is characterized 
analytically by the condition 
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(7) 



In Fig. y we depict by a solid line the ground-state en- 
ergy calculated according to (|^)-(0), for typical values of 
the parameters e = t' /t, t/ J and the hole concentration 
(1 — Ue)- We also depict by dashed lines the remnants 
of the individual crossing energy levels (|6|). Evidently, 
the ground-state energy (solid line) is an even periodic 
function of the external magnetic flux <&, with a macro- 
scopic energy barrier between different flux minima, in 
accordance with the Byers and Young |g^ characteriza- 
tion of a superconductor. The period is equal to $0 a-nd 
therefore the assignment q = 2e, discussed earlier on, 
leads to agreement with the observed flux quantization 
in the high-Tc superconducting copper oxide layers P6| . 
In order to establish flrmly the analytic result (H)-|3) 
we have also minimized numerically the classical energy 
(b)~®' Excellent agreement was obtained already for 
lattices with A = 20 x 20, as evidenced in the specific 
example of Fig. n^, where the open circles correspond to 
the numerical minimization data. 

Let us now turn our attention to the superfluid weight 
(or helicity modulus) Ds given by the curvature of the 
infinite lattice limit of the ground-state energy KE{^) at 
$ = [P5|J24|J2^, 



A, = A 



£0 

27r 



9$2 



Ds determines the ratio of the density of the superfluid 
charge carriers to their mass, and is related to the di- 
rectly measurable in-plane London penetration depth A^ 
by Ds = c2/(47re2A|). Quite generally, E{<^) has a 1/A^ 

expansion of the form £■($) = N"^ Eo{<^) + N Ei{<^) ^ 

which leads via (ph to a corresponding expansion for the 
superfluid weight Ds = N'^Df^ + ND^^^ H . Hence 



by exploiting the large- iV limit result (H)-(n) we get im- 
mediately the expression for the leading term Ds , 

Df^ ^U'n,{l-n,). (9) 

Our earlier arguments show that Ds is a measure of the 
stiffness of the classical phase fluctuations of the conden- 
sate. Furthermore, @j implies Ds = -Do, where Dq 
is the leading term in the 1/7V expansion of the Drude 
weight D = N^Do + NDi + ■■■ , studied in Ref. |l|] 
using the Kubo formalism for the current-current cor- 
relations. We have also verified, by a straightforward 
but lengthy calculation of £'i($) and the use of (g), that 
Ds = Di. Due to the analytic structure of the 1/A^ 
expansion, these results signify the term-by-term valid- 
ity of the identity Ds = D. Strictly speaking, of course, 
we have checked explicitly that Dg = D only up to and 
including terms Ds = Di, i.e., only up to and including 
leading quantum-fluctuation effects p8| . This, however, 
is sufficient for most practical purposes and permits us 
to exploit our calculations of the Drude weight, in the 
present study of the superfluid weight. For instance, the 
weight Ds ~ D, including leading quantum-fluctuation 
effects, is found [|l2[ to increase linearly with small hole 
concentration (1— ng) away from the half-fiUed-band limit 
{fie = 1). This trend, present already in (0), is a fun- 
damental characteristic of doped antiferromagnets. At 
higher doping values Dg = D eventually saturates and 
then starts to decrease. 

Concerning the expected transition temperature to the 
charged superfluid, i.e., superconducting, state under 
study we note that at a finite temperature T, the ratio of 
the thermal de Broglie wavelength of the charge carriers 
to their average distance is proportional to y^Ds/iksT), 
where Ds is the zero-temperature value determined by 
(^). Hence a naive application of the criterion for the 
occurrence of Bose-Einstein condensation in an ideal bo- 
son gas, whereby the latter ratio should become of order 
unity, suggests a transition temperature Tc of the form 



fcpT, = ADs 



(10) 



(8) 



$=0 



where A is a dimensionless constant of order unity. Of 
course, in the strictly two-dimensional model of continu- 
ous symmetry under study, a bona fide finite temperature 
phase transition can only be of the Kosterlitz-Thouless 
type which, nevertheless, leads again to an expression 
of the form dlQ). Indeed, the T/ii-structure of the clas- 
sical Hamiltonian (0)-(H) is a generalization of the two- 
dimensional XY model where the latter transition is well 



studied. In this context, it is important to note that a 
"universal" Hnear relation of the form ( pi3| ) has been es- 
tablished experimentally in the cuprates by Uemura et 
al. [pj in their remarkable study of Tc as a function of 
the zero-temperature value of A£ oc Dg. In the large- 
A'' limit, the Ds appearing in ( |l(] ) is just equal to Dg ' 
and the corresponding critical temperature Tc should 
be interpreted as the ordering temperature for the classi- 
cal phase fluctuations of the condensate, in analogy with 
the analysis of Emery and Kivelson l^^l of the classical 
phase fluctuations of the conventional BCS order param- 
eter. The higher order terms in the 1/A^ expansion of 
Ds — D capture the effects of the quantum fluctuations 
and renormalize downwards these weights [l2|, thereby 
reducing the corresponding value of Tc. 

Following the prescription of Emery and Kivelson [ p9[ , 
we have applied (HOh with A — 0.9; a numerical value 
extracted from the two-dimensional XY model 1 30 1 . Us- 
ing the calculated Dg ^ D oi Ref. |12], with the in- 
clusion of the leading quantum-fluctuation correction 
Ds = -Di , we depict in Fig. g the superconducting tran- 
sition temperature Tc as a function of the hole concentra- 
tion (1 — rig). Evidently, the dependence of Tc on (1 — rig) 
reflects that of Ds and reproduces qualitatively the ob- 
served trends in the cuprates pUTS]. In particular, we 
have that: Tc oc (1 — rie), for rig — s- 1. With an estimated 
J/ks ~ 1500 K in the cuprates |pl|, the value of Tc at 
optimum doping (1 — Uc) — 0.44 (0.36), seen in the solid 
(dashed) line of Fig. |, is Tc ~ 335 K (218 K). This pre- 
dicted value of Tc , signalling the coherence of the phase 
fluctuations of the condensate, should be regarded as an 
upper bound to an actual transition temperature because 
of the neglect of impurity disorder, higher-order quantum 
fluctuations, etc. From Fig. g we also note that with fur- 
ther hole doping Tc starts to decrease while beyond a 
critical doping value it vanishes, as the phase-modulated 
AF configuration, around which the present 1/N expan- 
sion is carried out, becomes unstable. 



IV. CONCLUSIONS 

We have demonstrated that flux quantization and a 
concomitant flnite value of superfluid weight Ds occur 
in the metallic phase-modulated AF ground state of the 
t-t'-J model (nh. By appealing to the universality class 
of the two-dimensional XY model, the corresponding su- 
perconducting transition temperature Tc is related to Ds 
linearly, via (IQ). The inclusion of leading quantum- 
fluctuation effects in Ds provides then a reasonable es- 
timate for the order of magnitude and the doping de- 
pendence of Tc in the cuprates. These results support 
our effective description of the charge carriers in terms 
of hard-core bosons. 
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FIG. 1. Ground-state energy vs external magnetic flux, 
for e = 0.45, t/J = 1.0 and 1 - n^ = 0.10. The zero-flux 
energy is subtracted to normalize the values. Solid line: the 
analytic result in the infinite lattice limit (A ~* oo), accord- 
ing to Eqs. (M)"(M). Dashed lines: remnants of the crossing 
energy-level parabolas discussed in the text. Open circles: nu- 
merical minimization results for the ground-state energy on 
a finite lattice (A = 20 x 20), as determined by Eqs. (|)-(0)- 
Evidently, the finite lattice numerical data (open circles) con- 
firm the infinite lattice limit analytic result (solid line). 



t/J=1.0 




FIG. 2. Superconducting transition temperature vs hole 
concentration, for t/J = 1.0 and e = 0.45 (solid line) or 
£ — 0.40 (dashed line), according to Eq. (Fui) with the in- 
clusion of leading quantum-fluctuation effects. 



